Variational cluster approach to correlated electron systems in low dimensions 
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A self-energy-functional approach is applied to construct cluster approximations for correlated 
lattice models. It turns out that the cluster-perturbation theory (Senechal et al, PRL 84, 522 
(2000)) and the cellular dynamical mean-field theory (Kotliar et al, PRL 87, 186401 (2001)) are 
limiting cases of a more general cluster method. Results for the one- dimensional Hubbard model 
are discussed with regard to boundary conditions, bath degrees of freedom and cluster size. 
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Low-dimensional systems of strongly interacting 
electrons, such as high-temperature superconductors, 
cuprate ladder compounds and organic conductors, cur- 
rently form a focus of intense experimental and theo- 
retical work. It is generally accepted that many of the 
fascinating physical properties of these materials arise 
from different kinds of short-range spatial correlations as 
well as from different phases with long-range order. On 
a low energy scale, this can be studied within effective 
models such as the Hubbard model [l|, for example. As 
these models have to be considered in the intermediate- 
to strong-coupling regime, weak-coupling perturbational 
approaches are inapplicable. Cluster methods, which ap- 
proximate the physics of the infinite system by solving 
the problem for a corresponding finite cluster, appear to 
be promising in this context as the interaction part can 
be treated numerically exact. Currently, there are two 
different groups of cluster methods which are discussed 
intensively: 

The point of reference for the first one is the direct clus- 
ter approach • Using an exact-diagonalization or quan- 
tum Monte-Carlo technique, the effect of short-range cor- 
relations can be studied by computing static and dynamic 
correlation functions for an isolated small cluster. This 
direct method, however, suffers from the fact that phase 
transitions and long-range order cannot occur in a sys- 
tem of finite size. Furthermore, the spectral function 
consists of a comparatively small number of poles. This 
has caused the recent development of an extension called 
cluster-perturbation theory (CPT) dHUH- Thc CPT 
procedure to calculate the one-electron Green's function 
G @ is sketched in Fig. for the Hubbard model: (i) 
The lattice is divided into small clusters, and the inter- 
cluster hopping V is switched off. The Green's function 
G' for the system of decoupled clusters (Hamiltonian H') 
is calculated numerically, (ii) The Green's function G of 
the lattice model H is then approximated by an RPA-like 
expression G = G'(l - VG')- 1 0. It has been pointed 
out that this corresponds to the first order in a systematic 
expansion in the inter-cluster hopping £| . 

The point of reference for the second type of clus- 
ter methods is the dynamical mean-field theory (DMFT) 



Thc lattice model H is mapped onto an impurity 
model H' consisting of a correlated site coupled to an infi- 
nite number of uncorrelated "bath" sites. The bath must 
be determined self-consistently. As a mean-field theory, 
the DMFT directly works in the thermodynamic limit 
and is, thus, able to describe phases with long-range or- 
der. Due to the locality of the self-energy however, it 
fails to incorporate the effects of short-range correlations. 
This has been the reason for the development of cluster 
extensions of the DMFT [UHfliEi. The main idea of 
the cellular DMFT (C-DMFT) Q is to replace the cor- 
related impurity site by a finite cluster (see Fig. Q and 
to proceed as follows: (i) The self-energy £ is calculated 
numerically for the system of decoupled clusters with an 
uncorrelated bath attached to each of the correlated sites 
(Hamiltonian H'). (ii) The (approximate) Green's func- 
tion of the lattice model H is then obtained from £ via 



the lattice Dyson equation: G = (G 



(iii) 



The parameters of the respective baths have to be recal- 
culated from G and X via the C-DMFT selfconsistency 
condition H[l3. This requires to repeat the above steps 
until self-consistency is reached. 

The first intention of this letter is to show that both 
the CPT and the C-DMFT can be considered as extreme 
limits of a more general cluster approach which is based 
on the self-energy-functional theory (SFT) proposed re- 
cently 15]. This answers an open question |5j for the 
relation between the different cluster methods and uni- 
fies two approaches which appear to be rather different 
at first sight. The presented cluster approach based on 
the SFT not only reproduces the CPT (n\, = 0) and 
the C-DMFT (rib = °°) but also allows to construct ap- 
proximations with an arbitrary number of bath sites rib- 
This intermediate approach does not waive any of the 
general merits such as causality and thermodynamical 
consistency. 

Secondly, this letter shows that there is room for 
new conceptual ideas beyond both the CPT and the C- 
DMFT: (i) A consistent cluster approach can be con- 
structed e.g. by attaching uncorrelated baths at the clus- 
ter boundary only, (ii) Not only the bath parameters but 
also the on-site energies of and thc hopping between the 
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system H 1 must share the same interaction part. For any 
X parameterized as S(t') we then have |l5| : 
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FIG. 1: Reference systems considered within various cluster 
approximations for the D = 2 Hubbard model. Filled cir- 
cles: on-site interaction. Solid lines: intra-cluster hopping. 
Dashed lines: inter-cluster hopping. Open circles: additional 
rib uncorrelated bath sites. Big circles: bath with nb = oo. 
Common to all methods is (i) the numerical solution of the 
system of decoupled clusters and (ii) the subsequent coupling 
of the clusters via an RPA-like or Dyson equation. Bath pa- 
rameters are determined self-consistently. The self-energy- 
functional theory (SFT) comprises the extreme limits CPT 
{ri h = 0) and C-DMFT (n b = oo). 



correlated sites may be determined in a self-consistent 
(variational) way. (iii) The question of boundary con- 
ditions can be answered by the method itself and need 
not be imposed by hand. To discuss the relevance and 
advantages or disadvantages of these points, some nu- 
merical results will be presented for the one-dimensional 
{D = 1) Hubbard model. 

Self-energy-functional theory (SFT): Consider a sys- 
tem of fermions on an infinite lattice with on-site 
Coulomb interactions at temperature T and chemical po- 
tential ii. Its Hamiltonian H = H^(t) + Hi(U) consists 
of a one-particle part which depends on a set of hop- 
ping parameters t and an interaction part with Coulomb- 
interaction parameters U. The grand potential fi can be 
obtained from the stationary point of a self-energy func- 
tional 



fi t [£] = TrlnHGo 1 - S)- 1 ) + F[£] 



(1) 



as has been discussed in Ref. [iff. Here Go = l/(u+fi— t) 
and F[S] is the Legendre transform of the Luttinger- 
Ward functional 3>[G]. As the latter is constructed as an 
infinite series of renormalized skeleton diagrams fl6| , the 
self-energy functional is not known explicitely. Neverthe- 
less, the exact evaluation of fit[X] an d the determination 
of the stationary point is possible jl5j on a restricted 
space S of trial self-energies S(t') € S. Due to this re- 
striction the procedure becomes approximative. 

Generally, the space S consists of t' representable self- 
energies. £ is termed t' representable if there are hopping 
parameters t' such that £ = £(i') is the exact self-energy 
of the model H' = H (t')+Hx(U) ("reference system"). 
Note that both the original system H and the reference 



fi t [£(t')] = r2' + Trln(-(Go 1 -S(t')) _1 )-Trhi( 



-G'), 
(2) 

where fi', G', and £(i') are the grand potential, the 
Green's function and the self-energy of the reference sys- 
tem H' while Go is the free Green's function of H. For a 
proper choice of t 1 (namely such that certain degrees of 
freedom in H' , e.g. those in different clusters, are decou- 
pled), a (numerically) exact computation of these quan- 
tities is possible. Hence, the self-energy functional @ 
can be evaluated exactly for this £ = £(i'). A certain 
approximation is characterized by a choice for £. As £ is 
parameterized by t' , this means to specify a space of vari- 
ational parameters t' . Any choice leads to a thermody- 
namically consistent approach since, once the variational 
procedure is carried out, Eq. (J2J provides an explicit ex- 
pression for a thcrmodynamical potential. For a further 
discussion of the general concept of the SFT see Ref. ; 
a detailed description of its practical application is given 
in Ref. [l^. 

Cluster approximations: Fig.^illustrates the construc- 
tion of cluster approximations within the framework of 
the SFT. To be specific, H is taken to be the Hubbard 
model with nearest-neighbor hopping. Subdividing the 
infinite lattice into identical clusters of finite size, H' is 
obtained from H by switching off the inter-cluster hop- 
ping and by switching on the hopping to new uncorre- 
lated (/7bath = 0) bath sites (Fig. [I] middle). Both oper- 
ations merely change the one-particle part of the Hamil- 
tonian, i.e. t — > t', while the interaction part (U) re- 
mains fixed - as required. To search for a stationary point 
on this space S of cluster-representable self-energies, one 
has to proceed as follows: (i) Compute the self-energy 
£(t') of the reference system for a given t' . (ii) Use 
Eq. © to evaluate fi t [£] at £ = £(*')• ( in ) Repeat 
steps (i) and (ii) for different t' to compute the function 
fi(i') = fit [£(£')] and the stationary point t' s given by 
dQ(t' s )/dt' = 0. As shown in Ref. [Fj|, causality require- 
ments are respected. 

The variational adjustment of the intra-cluster one- 
particle parameters t' can be looked upon as a (partial) 
compensation for the error introduced by the finite clus- 
ter size. An inclusion of rib bath sites per original cor- 
related site enlarges the number of variational parame- 
ters and thereby the space S. This is expected to (and 
does) improve the approximation (see results below). In 
the limit of infinite cluster size (number of correlated 
sites within a cluster N c — > oo), the exact self-energy 
becomes t' representable and therefore the cluster ap- 
proximation itself becomes exact. Since N c must be fi- 
nite (small) in any practical calculation, one should focus 
on local quantities such as the on-site Green's function 
G tl = (Gq 1 - £(*£))«S for example. 

CPT and C-DMFT: For a certain cluster approxima- 
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tion, it has to be specified which of the different intra- 
cluster one-particle parameters t' are treated as vari- 
ational parameters. The simplest idea is to consider 
the intra-cluster hopping as fixed at the original values, 
iL = tij (for i, j in the same cluster), and not to switch 
on a hopping to bath sites (i. e. not to introduce any bath 
sites). In this case there is no variational parameter at 
all. E(t') is calculated once, and the Green's function G 
for the original model is obtained by the (lattice) Dyson 
equation G 1 = G^ 1 - S(t'). As £(*') is the exact self- 
energy for H' we also have: S(t') = Gq _1 — G' _1 and 
consequently G 1 = G^ 1 - G'^ 1 + G'' 1 = G 1 ' 1 - V 
with V being the inter-cluster hopping. As this is equiv- 
alent to the RPA-type equation mentioned above, one 
recovers the CPT. 

The C-DMFT is obtained by introducing a hopping 
t' ir to rib = oo bath sites r = l,...,rib per correlated 
site i = 1,...,N C and taking this hopping ("hybridiza- 
tion" ) and the bath on-site energies e£ r as variational 
parameters while for the correlated sites t'^ = tij is still 
fixed. Assume that bath parameters {t' ir ,e' ir } can be 
found such that the C-DMFT self-consistency equation 
is fulfilled. In Ref. this was given in fc-space represen- 
tation (k from the reduced Brillouin zone). In the real- 
space representation the self-consistency equation reads 
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(G 1 — S({^>, e'i r }))i j — , 



(3) 



where i,j must belong to the same cluster. This imme- 
diately implies that H({t' ir , e' ir }) satisfies the SFT Euler 
equation d£l(t')/dt' = or, calculating the derivative, 



i 



Gq - S(i') 



- G' 



dt' 



= 0. 



(4) 



This holds since the "projector" dHij(t') / 'dt' = if i,j 
belong to different clusters as these are decoupled in the 
reference system. We conclude that the self-energy func- 
tional is stationary at the C-DMFT self-energy. 

In principle, approximations may also be constructed 
in reciprocal k space. For the Hubbard model, however, 
there is no simple reference system as the interaction part 
is non-local in k space. In particular, it is not possible to 
recover the dynamical cluster approximation (DCA) 
within the SFT. 

"Intermediate" approach: The C-DMFT self- 
consistency equation can generally be fulfilled for 
?ib = oo only. Within a cluster approach based on 
the SFT there arc no formal problems, however, if 
rib < oo. A finite yields an approximation inferior 
as compared to rib = oo (C-DMFT) and superior as 
compared to rib = (CPT) as there are less or more 
variational parameters, respectively. For < rib < oo 
the parameters must be found to satisfy the Eq. Q. 

Since the convergence with respect to rib appears to 
be rapid for local physical quantities [lij and since the 
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FIG. 2: Eo = Q + fJ-(N) vs. t' as obtained by evaluating 
the self-energy functional Q — fi[E(t')]. Original system H: 
D = 1 Hubbard model with n.n. hopping t = 1 and U = 8 
for T — and y, = {7/2 (half-filling). Reference system H': 
set of decoupled Hubbard chains with N c sites each; N c = 
2, 4, 6, 8, 10 as indicated. Variational parameter: n.n. hopping 
t' of H' . The inset shows AEq — Eo — So, mm vs. t' . 



cluster Hilbert-space dimension increases exponentially 
with rib, approximations with small ny, (or even rib = 0, 
CPT) appear advantageous. Anyway, the bath concept 
must become irrelevant in the limit N c —> oo. On the 
other hand, there are good reasons to introduce bath 
sites: Depending on the dimensionality of the problem, 
it can be the description of the local (temporal) correla- 
tions that needs to be improved in first place. Note that 
for N c — 1 and rib = oo the DMFT is recovered which 
represents the exact solution for D = oo 0- Further- 
more, bath sites can serve as particle reservoirs which 
will be essential for a proper description of filling de- 
pendencies. Finally, the presence of bath sites may also 
facilitate practical calculations to treat the reference sys- 
tem, e.g. by an attenuation of the sign problem in the 
context of a Hirsch-Fye like QMC approach . 

Intra-cluster hopping determined variationally: For 
both the CPT and the C-DMFT the hopping between 
correlated sites is fixed at ty = tij (for the C-DMFT, this 
is even a necessary condition to satisfy Eq. as can eas- 
ily be seen from a high-frequency expansion). Contrary, 
within the SFT there is a priori no reason to fix t\j. A 
cluster approximation with £L determined variationally 
represents another "intermediate" approach as shall be 
discussed in the following. 

Numerical calculations have been performed for the 
D = 1 Hubbard model. Instead of solving the Euler 
equation we have directly evaluated the self-energy 
functional according to Eq. (J5J. The reference system H' 
is taken to be a set of decoupled Hubbard chains with 
A^c correlated sites each. For 7V C < 10 (no additional 
bath sites), the ground-state energy E' Q (and thereby 
Q' = E' - fj.(N)' for T = 0) and the Green's function G' 
are computed using the standard Lanczos algorithm |l9| . 



The self-energy is obtained as E(t') = G 



G' 



In 



4 



-0.20 
-0.22 
-0.24 
-0.26 
-0.28 
-0.30 
-0.32 
-0.34 



1 1 1 1 1 




1 1 1 1 1 1 1 
*r / 






C \ \- 

/ 


t t t r 




t t t r 


N =2/ " 

c / 

/4/ 




N =4/ / 

L C / / 

V 6J/7 

v\ 8 


/ / Q. 




\ \ -in JrJf 

\\\ ^"-j. n - 














1,1,1 




1,1,1,1 



direct: 
i-N =2 

c 



-6 

- exact 



0.2 0.4 0.6 -1 -0.5 0.5 1 



t 



t 



FIG. 3: Eo as in Fig. |5]but vs. different variational parame- 
ters: n.n. hopping th to uncorrelated bath sites (left), hopping 
t r between the edge sites (right). Dashed lines: Eo without 
parameter optimization. Arrows: Eo calculated for an iso- 
lated cluster. Exact result for Eo from Ref. |2CB . 



Eq. (J2J , the trace "Tr" consists of a sum over a complete 
set of one-particle quantum numbers and (after analytical 
continuation to the real u axis) a frequency integration. 
To keep the calculations simple, only a single variational 
parameter is taken into account. 

Fig. shows E (t') = Q(t') + fi(N)' with fi(t') = 
f2t[S(i')] where t' is the nearest-neighbor hopping within 
the cluster. Clearly, Eq is stationary for t' very close to 
t but t' ^ t. The effect is most obvious for the small- 
est cluster size N c — 2 (see inset). We conclude that 
a variational determination of the hopping between the 
correlated sites in fact improves the approximation. Sur- 
prisingly, however, the improvement is almost negligible 
for reasonable N c . Similar results are obtained when dif- 
ferent selected hopping parameters are varied. 

Different variational parameters: It is possible to con- 
struct a consistent cluster approach by attaching uncorre- 
lated bath sites only at the boundaries of the respective 
cluster. Compared to the C-DMFT, this choice repre- 
sents a certain restriction of the parameter set which is 
motivated by the expectation that here bath degrees of 
freedom compensate for the finite-size errors most effi- 
ciently. Fig. |21 (left) shows that switching on the hopping 
tb to two bath sites in fact lowers the minimal Eq. With 
increasing N c the energy difference E (tb = t™ in ) — E o (0) 
decreases, and E (tb) becomes almost flat for small tb, 
as expected. Note that the binding energy gain due to 
inclusion of two bath sites \E (N C , t^ in ) - E o (N c ,0)\ is 
always smaller than the gain \E (N C + 2, 0) - E (N C , 0)| 
due to a larger cluster: Introducing bath sites is less ef- 
ficient as increasing the cluster size, at least for D = 1. 
Interestingly, the convergence of Eq with increasing 7V C 
appears to be much faster as compared to the direct clus- 
ter method (see arrows in Fig.|3). 

Within the CPT the question of the correct boundary 



conditions is decided a posteriori by inspection of the 
respective results for the spectral density [f|. Here we 
introduce a hopping parameter i r between the edge sites 
and let the method "decide" by itself. As is seen in Fig. 
13 (right), a minimum for Eq is obtained at t T — (open 
boundary conditions) while there is no indication for a 
stationary point at t T = t (periodic boundary conditions) . 

Conclusion: The self-energy-functional theory has 
been shown to provide a consistent and very general 
framework for the construction of variational cluster ap- 
proximations including the CPT and the C-DMFT. The 
extreme flexibility with regard to the choice of variational 
parameters offers a variety of further methodical devel- 
opments and applications. 
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